tions  reveals  well-preserved  plant  organs 
and  tissues,  with  conifer  wood  possess¬ 
ing  distinctive  secondary  xylem  predom¬ 
inating.  Young  roots  with  well-devel¬ 
oped  tracheids  and  wood  rays  reveal  a 
tissue  surrounding  the  xylem,  which 
may  be  phloem  with  definite  alternating 
rows  of  thick-walled  fibers  (Fig.  IE). 
The  cortex  is  usually  not  well  preserved 
(Fig.  1C),  whereas  the  epidermis  when 
present  is  sharply  delineated  by  its 
rather  wavy  cuticular  layer  (Fig.  1C). 
The  tracheids  of  the  primary  xylem 
possess  a  helical  secondary  wall;  those 
of  the  secondary  xylem  have  uniseriate, 
bordered  pits  with  crassulae  on  the 
radial  section  (Fig.  1G).  The  rays  are 
uniseriate,  ranging  from  two  to  eight 
cells  (Fig.  II). 

Other  structures  observed  include  a 


seed  of  possible  conifer  origin  with  a 
well-preserved  embryo  and  several  well- 
preserved  fern  annuli. 
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Around-the-World  Atomic  Clocks: 

Predicted  Relativistic  Time  Gains 

Abstract.  During  October  1971,  four  cesium  beam  atomic  clocks  were  flown 
on  regularly  scheduled  commercial  jet  flights  around  the  world  twice ,  once  east¬ 
ward  and  once  westward,  to  test  Einstein's  theory  of  relativity  with  macroscopic 
clocks.  From  the  actual  flight  paths  of  each  trip,  the  theory  predicts  that  the 
flying  clocks,  compared  with  reference  clocks  at  the  U.S.  Naval  Observatory, 
should  have  lost  40  ±  23  nanoseconds  during  the  eastward  trip,  and  should  have 
gained  275  ±  21  nanoseconds  during  the  westward  trip.  The  observed  time  differ¬ 
ences  are  presented  in  the  report  that  follows  this  one. 


One  of  the  most  enduring  scientific 
debates  of  this  century  is  the  relativistic 
clock  “paradox”  (2)  or  problem  (2), 
which  stemmed  originally  from  an 
alleged  logical  inconsistency  in  predicted 
time  differences  between  traveling  and 
reference  clocks  after  a  round  trip. 
This  seemingly  endless  theoretical  de¬ 
bate,  which  has  flared  up  recently  with 
renewed  vigor  (2,  3),  begs  for  a  con¬ 
vincing  empirical  resolution  with  macro¬ 
scopic  clocks.  A  simple  and  direct  ex¬ 
perimental  test  of  the  clock  problem 
with  portable  atomic  clocks  is  now  pos¬ 
sible  because  of  the  unprecedented 
stability  achieved  with  these  clocks  ( 4 ). 

In  this  first  of  two  reports,  we  pre¬ 
sent  relativistic  time  differences  calcu¬ 
lated  from  flight  data  for  our  recent 
around-the-world  flying  clock  experi¬ 
ments.  The  theory  predicts  a  detectable 
effect  with  cesium  beam  clocks  if  they 
are  flown  around  the  world  at  typical 
jet  aircraft  speeds  ( 4 ).  Moreover,  it 
predicts  an  interesting  asymmetry  in  the 
time  difference  between  the  flying 
clocks  and  a  ground  reference  clock, 
depending  on  the  direction  of  the  cir¬ 
cumnavigation  ( 4 ).  Predicted  time  dif¬ 


ferences  are  compared  with  our  ob¬ 
served  time  differences  in  the  following 
report. 

A  brief  elementary  review  of  the 
theory  seems  appropriate,  particularly 
because  of  some  confusion  about  the 
capacity  of  such  experiments  to  pro¬ 
duce  meaningful  results  (5).  Special 
relativity  predicts  that  a  moving  stan¬ 
dard  clock  will  record  less  time  com¬ 
pared  with  (real  or  hypothetical)  co¬ 
ordinate  clocks  distributed  at  rest  in 
an  inertial  reference  space.  For  low 
coordinate  speeds  (u2  <  c2),  the  ratio 
of  times  recorded  by  the  moving  and 
reference  coordinate  clocks  reduces  to 
(1  —  u2/2c2),  where  c  is  the  speed  of 
light.  Because  the  earth  rotates,  stan¬ 
dard  clocks  distributed  at  rest  on  the 


Table  1.  Predicted  relativistic  time  differences 
(nsec). 


Effect 

Direction 

East 

West 

Gravitational 

144  ±  14 

179  ±  18 

Kinematic 

—184  ±  18 

96  ±  10 

Net 

-40  ±  23 

275  ±  21 

surface  are  not  suitable  in  this  case  as 
candidates  for  coordinate  clocks  of  an 
inertial  space.  Nevertheless,  the  relative 
timekeeping  behavior  of  terrestrial  clocks 
can  be  evaluated  by  reference  to  hypo¬ 
thetical  coordinate  clocks  of  an  under¬ 
lying  nonrotating  (inertial)  space  (6). 

For  this  purpose,  consider  a  view  of 
the  (rotating)  earth  as  it  would  be 
perceived  by  an  inertial  observer  look¬ 
ing  down  on  the  North  Pole  from  a 
great  distance.  A  clock  that  is  station¬ 
ary  on  the  surface  at  the  equator  has  a 
speed  RQ  relative  to  nonrotating  space, 
and  hence  runs  slow  relative  to  hypo¬ 
thetical  coordinate  clocks  of  this  space 
in  the  ratio  1  —  R2Q2/2c2,  where  R 
is  the  earth’s  radius  and  Q  its  angular 
speed.  On  the  other  hand,  a  flying  clock 
circumnavigating  the  earth  near  the 
surface  in  the  equatorial  plane  with  a 
ground  speed  v  has  a  coordinate  speed 
RQ  +  v,  and  hence  runs  slow  with  a 
corresponding  time  ratio  1  —  (RQ  4- 
v)2/2 c2.  Therefore,  if  r  and  r0  are  the 
respective  times  recorded  by  the  flying 
and  ground  reference  clocks  during  a 
complete  circumnavigation,  their  time 
difference,  to  a  first  approximation,  is 
given  by 


r  -  to  =  -  (2 RQv  +  v2)  ro/2c2  ( 1 ) 

Consequently,  a  circumnavigation  in  the 
direction  of  the  earth’s  rotation  (east¬ 
ward,  i;  >  0)  should  produce  a  time 
loss,  while  one  against  the  earth’s  rota¬ 
tion  (westward,  v  <  0)  should  produce  a 
time  gain  for  the  flying  clock  if  ]v|  - — ■ 
RQ. 

General  relativity  predicts  another 
effect  that  (for  weak  gravitational  fields) 
is  proportional  to  the  difference  in  the 
gravitational  potential  for  the  flying  and 
ground  reference  clocks.  If  the  surface 
value  of  the  acceleration  of  gravity  is  g 
and  the  altitude  for  the  circumnaviga¬ 
tion  is  h  <  R,  the  potential  difference  is 
gh,  and  Eq.  1  then  reads 


r—ro  =  Igh/C2 


(2Rttv  +  V2)/2c2]r0 

(2) 


The  gh/c 2  term,  which  is  related  to  the 
gravitational  “red  shift,”  predicts  a  time 
gain  for  the  flying  clock  irrespective  of 
the  direction  of  the  circumnavigation. 
For  typical  aircraft  speeds  and  altitudes, 
both  the  gravitational  and  kinematic 
terms  in  Eq.  2  are  comparable  in  abso¬ 
lute  magnitude,  and  v2/  2c2  is  small 
compared  with  RQv/c2.  For  a  westward 
circumnavigation  (v  <  0)  both  terms 
are  positive  and  they  add  to  give  a  large 
net  time  gain,  but  for  an  eastward  cir¬ 
cumnavigation  (v  >  0)  they  tend  to 
cancel  and  produce  a  net  time  differ- 
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Fig.  1.  Predicted  relativis¬ 
tic  time  gain  for  a  flying 
clock  after  a  nonstop 
equatorial  circumnaviga¬ 
tion  of  the  earth  at  var¬ 
ious  altitudes.  The  label¬ 
ed  dots  correspond  to 
cruising  ground  speeds 
and  altitudes  for  the  in¬ 
dicated  aircraft.  The  area 
within  the  hatched  lines 
is  below  approximate  de¬ 
tection  thresholds  with  a 
portable  cesium  beam 
clock. 


ence  that  may  be  positive  or  negative 
depending  on  details  of  the  flight. 

We  can  compare  the  predicted  time 
differences  with  detection  thresholds.  If 
the  circumnavigation  is  nonstop,  to  a 
first  approximation  the  trip  time  t0  = 
27tJR/|v|.  Substitution  of  this  value  for 
,t0  in  Eq.  2  gives 

r  -  r„=^  [gh/\v\  -  IHlv/|v|  -  |v|/2] 

(3) 

This  relationship  is  graphically  illus¬ 
trated  in  Fig.  1  over  the  range  of 
ground  speeds  and  altitudes  of  interest. 
The  area  within  the  hatched  lines  in 
Fig.  1  is  below  detection  thresholds 
estimated  from  past  experience  at  the 
U.S.  Naval  Observatory  with  portable 
cesium  beam  clocks  (7).  The  labeled 
points  correspond  to  cruising  altitudes 
and  ground  speeds  for  the  indicated 
aircraft  ( 8 ).  Figure  1  illustrates  that 
circumnavigations  with  cesium  beam 
clocks  at  jet  aircraft  speeds  should  pro¬ 
duce  measurable  relativistic  time  dif¬ 
ferences.  Furthermore,  the  mere  ex¬ 
istence  of  a  definite  east-west  directional 
asymmetry  in  the  observed  time  dif¬ 
ferences  would  give  strong  evidence  for 
the  validity  of  the  kinetic  term  in  Eq.  2. 

These  predictions  for  equatorial  cir¬ 
cumnavigations  will  be  modified  some¬ 
what  for  our  actual  flights  because 
commercial  around-the-world  jet  flights, 
of  course,  do  not  follow  equatorial 
paths,  nor  do  they  maintain  constant 
altitude,  ground  speed,  or  latitude.  In 
this  case,  however,  it  is  only  necessary 
to  integrate  the  appropriate  differential 
form  of  Eq.  2  along  the  actual  flight 
path: 

r  —  To  — 

[ghf c~  —  (IRilvcosdcosX  -f  v2)/2c2]  dr 

(4) 

This  expression  contains  a  slightly  mod¬ 
ified  directionally  dependent  term, 
which  for  nonequatorial  flights  becomes 
proportional  to  both  the  eastward  com¬ 
ponent  of  the  ground  speed,  v  •  cos  0, 
and  the  cosine  of  the  latitude,  cos  X 
(4,  9).  Because  only  lowest  order  rela¬ 
tivistic  time  differences  can  be  detected 
in  our  experiments,  only  lowest  order 
terms  need  be  retained  in  the  calculated 
predictions,  and,  to  this  order  of  ap¬ 
proximation,  it  is  immaterial  whether 
the  differential  time  for  r  or  for  r0  is 
used  in  Eq.  4. 

The  eastward  trip  began  on  4  October 
1971  at  19h30m  U.T.  and  lasted  65.4 
hours  with  41.2  hours  in  flight.  The 
westward  trip  began  during  the  follow¬ 
ing  week  on  13  October  at  19h40m  U.T. 


and  lasted  80.3  hours  with  48.6  hours  in 
flight.  Flight  data  necessary  for  numeri¬ 
cal  evaluation  of  Eq.  4  were  provided 
by  the  various  flight  captains.  In  most 
cases  they  traced  their  flight  path  on 
an  appropriate  flight  map  and  recorded 
the  time  and  aircraft  ground  speed  and 
altitude  at  various  navigation  check 
points  along  the  flight  path.  This  in¬ 
formation  divided  the  eastward  trip  into 
125  intervals  and  the  westward  trip  into 
108  intervals.  The  latitude  and  longi¬ 
tude  for  each  check  point,  read  di¬ 
rectly  from  the  flight  maps,  combined 
with  the  time  (U.T.)  over  each  check 
point  permits  calculation  of  an  average 
ground  speed,  latitude,  and  eastward 
azimuth  for  each  interval.  The  average 
altitude  for  each  interval  was  taken  as 
the  average  of  the  altitudes  at  the  end 
points.  This  information  then  permits 
numerical  evaluation  of  the  integral  in 
Eq.  4.  Table  1  gives  the  predicted  time 
differences  resulting  from  these  calcu¬ 
lations. 

We  conclude  this  report  with  a  word 
about  uncertainty  in  these  predictions. 
Possible  errors  stem  from  two  sources: 
(i)  errors  and  deficiencies  in  the  flight 
data,  and  (ii)  theoretical  approximations 
used  in  the  derivation  of  Eq.  4.  We  es¬ 
timate  the  maximum  possible  fractional 
uncertainty  from  flight  data  errors  to 
be  less  than  1 0  percent  for  each  term  of 
Eq.  4  after  numerical  integration.  This 
estimate  includes  both  systematic  and 
random  errors  in  the  flight  data.  If  it 
is  assumed  that  the  errors  from  these 
terms  add  in  quadrature,  the  maximum 
fractional  uncertainty  in  the  net  value 
for  the  eastward  circumnavigation  is 
about  60  percent,  while  that  for  the 
westward  circumnavigation  is  only  8 
percent.  These  uncertainties  are  listed  in 
Table  1. 

Although  neglect  of  higher  order 
terms  (proportional  to  c~4,  c~f>,  and 
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so  forth)  in  our  theoretical  approxima¬ 
tions  is  fully  justified,  small  but  perhaps 
not  entirely  negligible  first  order  effects 
may  arise  from  the  presence  of  the 
moon  and  sun.  In  fact,  the  center-of- 
mass  of  the  earth-moon  system,  not  the 
center  of  the  earth,  is  in  free  fall 
around  the  sun,  and  a  more  precise 
calculation  should  include  this  effect.  It 
is  unlikely,  however,  that  the  precision 
of  our  experiments  permits  detection  of 
any  effects  other  than  the  dominating 
ones  retained  in  Eq.  4  ( 10 ,  11). 
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Around-the-World  Atomic  Clocks: 
Observed  Relativistic  Time  Gains 


Abstract.  Four  cesium  beam  clocks  flown  around  the  world  on  commercial 
jet  flights  during  October  1971,  once  eastward  and  once  westward ,  recorded 
directionally  dependent  time  differences  which  are  in  good  agreement  with  predic¬ 
tions  of  conventional  relativity  theory.  Relative  to  the  atomic  time  scale  of  the 
U.S .  Naval  Observatory,  the  flying  clocks  lost  59  ±10  nanoseconds  during  the 
eastward  trip  and  gained  273  ±7  nanoseconds  during  the  westward  trip ,  where  the 
errors  are  the  corresponding  standard  deviations.  These  results  provide  an  un¬ 
ambiguous  empirical  resolution  of  the  famous  clock  “ paradox' 9  with  macroscopic 
clocks. 


In  science,  relevant  experimental  facts 
supersede  theoretical  arguments.  In  an 
attempt  to  throw  some  empirical  light 
on  the  question  of  whether  macroscopic 
clocks  record  time  in  accordance  with 
the  conventional  interpretation  of  Ein¬ 
stein’s  relativity  theory  (J),  we  flew 
four  cesium  beam  atomic  clocks  around 
the  world  on  commercial  jet  flights, 
first  eastward,  then  westward.  Then  we 
compared  the  time  they  recorded  dur¬ 
ing  each  trip  with  the  corresponding 
time  recorded  by  the  reference  atomic 
time  scale  at  the  U.S.  Naval  Observa¬ 
tory,  MEAN(USNO)  (2).  As  was  ex¬ 
pected  from  theoretical  predictions  (I), 
the  flying  clocks  lost  time  (aged  slower) 
during  the  eastward  trip  and  gained 
time  (aged  faster)  during  the  westward 
trip.  Furthermore,  the  magnitudes  of 
the  time  differences  agree  reasonably 
well  with  predicted  values,  which  were 
discussed  in  the  preceding  report  ( 1 ) .  In 
this  second  report,  we  present  the  time 
difference  data  for  the  flying  ensemble, 
and  explain  the  methods  by  which  the 
relativistic  time  differences  were  ex¬ 
tracted. 

The  development  of  compact  and 
portable  cesium  beam  atomic  clocks  ( 3 ) 
permits  a  terrestrial  test  of  relativity 


theory  with  flying  clocks.  The  funda¬ 
mental  unit  of  time  interval,  the  sec¬ 
ond,  is  now  by  definition  equal  to 
9,192,631,770  accumulated  periods  of 
the  frequency  of  the  atomic  transitions 
of  an  “ideal”  cesium  beam  frequency 
standard  (2,  3).  Because  these  clocks 
are  regulated  by  the  frequency  of  a 
natural  atomic  transition,  a  particularly 
well  defined  hyperfine  transition  in  the 
ground  state  of  the  133Cs  atom,  they 

Table  1.  Observed  relativistic  time  differences 
from  application  of  the  correlated  rate-change 
method  to  the  time  intercomparison  data  for 
the  flying  ensemble.  Predicted  values  are  listed 
for  comparison  with  the  mean  of  the  ob¬ 
served  values;  S.D.,  standard  deviation. 


Clock 

serial 

No. 

At  (nsec) 

Eastward* 

Westward 

120 

-  57 

277 

361 

-74 

284 

408 

-  55 

266 

447 

-  51 

266 

Mean 

±  S.D. 

-  59  ±  10 

273  ±7 

Predicted 

±  Error  est. 

-  40  ±  23 

275  ±  21 

*  Negative  signs  indicate  that  upon  return  the 
time  indicated  on  the  flying  clocks  was  less  than 
the  time  indicated  on  the  MEAN(USNO)  clock 
of  the  U.S.  Naval  Observatory. 


approach  the  ideal  standard  clock  of 
relativity  theory. 

However,  no  two  “real”  cesium  beam 
clocks  keep  precisely  the  same  time, 
even  when  located  together  in  the  labo¬ 
ratory,  but  generally  show  systematic 
rate  (or  frequency)  differences  which  in 
extreme  cases  may  amount  to  time  dif¬ 
ferences  as  large  as  1  /x sec  per  day.  Be¬ 
cause  the  relativistic  time  offsets  expect¬ 
ed  in  our  experiments  are  only  of  the 
order  of  0.1  jtsec  per  day  (1,  4),  any 
such  time  divergences  (or  rate  differ¬ 
ences)  must  be  taken  into  account. 

A  much  more  serious  complication  is 
caused  by  the  fact  that  the  relative  rates 
for  cesium  beam  clocks  do  not  remain 
precisely  constant.  In  addition  to  short 
term  fluctuations  in  rate  caused  mainly 
by  shot  noise  in  the  beam  tubes,  cesium 
beam  clocks  exhibit  small  but  more 
or  less  well  defined  quasi-permanent 
changes  in  rate.  The  times  at  which 
these  rate  changes  occur  typically  are 
separated  by  at  least  2  or  3  days  for 
good  clocks.  Some  clocks  have  been  ob¬ 
served  in  the  laboratory  to  go  as  long  as 
several  months  without  a  rate  change 
(2,  5). 

These  unpredictable  changes  in  rate 
produce  the  major  uncertainty  in  our 
results.  Because  of  the  nature  of  these 
changes,  however,  their  effect  on  the 
observed  time  differences  can  be  re¬ 
moved  to  a  large  extent  in  the  data 
analysis.  Under  normal  conditions 
changes  in  relative  rates  occur  inde¬ 
pendently,  that  is,  there  are  no  known 
systematic  correlations  between  rate 
changes  of  one  clock  and  those  of  an¬ 
other.  Consequently,  the  chance  that 
two  or  more  clocks  will  change  rate  by 
the  same  amount  in  the  same  direction 
at  the  same  time  is  extremely  remote. 
Because  of  the  random  and  independent 
character  of  these  rate  changes,  the 
long-term  average  rate  of  an  ensemble 
of  clocks  is  more  stable  than  the  rate 
of  any  individual  member. 

Starting  at  0h  U.T.  on  25  September 
1971,  we  recorded  more  than  5000  time 
differences  during  the  data  period. 
Figure  1  shows  the  time  difference  data 
relative  to  MEAN(USNO)  for  the  en¬ 
tire  data  period,  which  lasted  636  hours. 
The  labels  in  Fig.  1  are  the  serial  num¬ 
bers  of  the  corresponding  clocks,  and 
the  traces  give  the  measured  differences 
in  time  between  the  corresponding 
clock  and  MEAN(USNO).  Of  course, 
no  comparisons  with  MEAN(USNO) 
were  possible  during  the  trips.  Exactly 
the  same  electronic  arrangement  was 
used  for  all  time  intercomparison  mea- 
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surements,  and  time  differences  were 
measured  with  an  electronic  time  inter¬ 
val  counter  to  the  nearest  nanosecond 
(5). 

The  axis  of  abscissas  in  Fig.  1  is  the 
running  time  for  the  data  period,  and, 
because  at  any  given  instant  all  clocks 
involved  read  the  same  time  to  within 
much  better  than  the  nearest  second, 
it  is  immaterial  whether  the  running 
time  is  the  time  of  MEAN(USNO)  or 
any  member  of  the  flying  ensemble.  The 
times  at  which  the  time  differences  were 
measured  are  certain  to  within  about 
±  5  minutes. 

The  tangential  slope  of  a  trace  in 
Fig.  1  equals  the  instantaneous  relative 
rate  for  the  corresponding  clock.  For 
example,  a  positive  slope  of  30°  means 
that  the  clock  was  gaining  on  MEAN- 
(USNO)  by  about  10  nsec  per  hour. 
Figure  1  shows  that  no  member  of  the 
flying  ensemble  maintained  a  constant 
rate  relative  to  MEAN(USNO)  during 
the  entire  data  period.  Spontaneous  rate 
changes  between  periods  of  virtually 
constant  rate  are  apparent. 

Figure  2  gives  magnified  views  of  the 
average  time  difference  data  of  Fig.  1 
immediately  before  and  after  each  trip. 
These  magnified  views  clearly  show 
that,  on  the  average,  the  flying  clocks 
lost  time  during  the  eastward  trip  and 
gained  time  during  the  westward  trip. 
Hence  the  predicted  east-west  direc¬ 
tional  asymmetry  (1)  is  qualitatively 
confirmed. 

We  have  extracted  relativistic  time 
differences  from  our  data  with  two  dif¬ 
ferent  methods,  the  “average  rate” 
method  and  the  “correlated  rate-change” 
method  (6).  Although  the  average  rate 
method  is  less  reliable,  it  is  the  simpler 
of  the  two  methods  and  we  will  describe 
it  first. 

Figure  2  shows  least-squares  straight- 
line  fits  to  the  average  data  of  Fig.  1 
for  an  interval  of  25  hours  immediately 
before  and  after  each  trip.  The  slopes 
for  these  lines  give  the  mean  rates  for 
this  interval  (electronic  noise  causes  the 
short-term  fluctuations  about  the  local 
mean  rate).  If  we  assume  that  only 
one  rate  change  occurred  in  the  “aver¬ 
age”  clock  during  each  trip,  and  that 
it  occurred  at  the  midpoint,  relativistic 
time  differences  follow  from  extrapola¬ 
tion  with  the  average  of  the  initial  and 
final  rates.  The  slopes  of  the  dashed 
extrapolation  lines  in  Fig.  2  are  equal 
to  the  average  of  the  slopes  of  the  linear 
fits,  and  with  this  fitting  interval  the 
average  rate  method  gives  a  time  loss 
of  66  nsec  for  the  eastward  trip  and  a 


Fig.  1.  Time  differences 
between  each  clock  in  the 
flying  ensemble  and  the 
Naval  Observatory  clock 
MEAN(USNO)  at  hourly 
intervals  throughout  the 
data  period.  The  hour 
count  along  the  abscissa 
begins  at  0h  U.T.,  25  Sep¬ 
tember  1971.  The  trace  of 
each  clock  is  labeled  with 
the  clock’s  serial  number; 
and  the  trace  labeled  Av¬ 
erage  is  the  average  of 
the  four  time  differences. 


time  gain  of  205  nsec  for  the  westward 
trip. 

Reliability  of  results  with  the  average 
rate  method,  however,  depends  on  the 
unlikely  chance  that  only  one  rate 
change  occurred  during  each  trip  and 
that  it  occurred  at  the  midpoints. 
Furthermore,  there  is  no  obvious 
method  for  estimating  the  experimental 
error.  Nevertheless,  the  average  rate 
method  does  produce  convincing  quali¬ 
tative  results. 

The  correlated  rate-change  method 
produces  more  reliable  quantitative  re¬ 
sults  because,  with  this  method,  rate 
changes  for  individual  members  of  the 
flying  ensemble  during  each  trip  are 
taken  into  account.  Because  rate  changes 
relative  to  MEAN(USNO)  are  uncorre¬ 
lated  whether  or  not  there  is  an  actual 
comparison  with  MEAN(USNO),  rec¬ 
ords  of  the  time  differences  between 


each  member  of  the  flying  ensemble  per¬ 
mit  reconstruction  of  the  rate-change 
history  for  each  clock  during  the  trips. 

The  correlated  rate-change  method  is 
basically  an  old  technique  (7).  It  was 
first  used  with  atomic  clocks  to  supple¬ 
ment  a  frequency  averaging  method 
that  employed  very  low  frequency  radio 
transmissions  (5),  and  then  became  a 
main  element  in  the  generation  of  the 
Naval  Observatory  atomic  time  scale. 
The  computer  algorithm  (2)  currently 
used  to  construct  MEAN(USNO) 
evolved  from  studies  of  time  scales  con¬ 
structed  with  the  correlated  rate-change 
method.  Because  the  stability  of  a  single 
cesium  beam  clock  would  probably  be 
inadequate  to  permit  an  absolutely  un¬ 
ambiguous  detection  of  the  expected 
relativistic  effects,  we  used  four  flying 
clocks  so  that  the  correlated  rate-change 
or  other  similar  methods  could  be  ap- 
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Fig.  2.  Magnified  views  of  the  average  time  difference  data  of  Fig.  1  in  the  vicinity 
of  (a)  the  eastward  trip  and  (b)  the  westward  trip.  The  small  step  down  for  the  east¬ 
ward  trip  indicates  an  average  time  loss,  and  the  somewhat  larger  step  up  for  the 
westward  trip  indicates  an  average  time  gain.  The  indicated  time  differences  follow 
from  application  of  the  average  rate  method  with  a  fitting  interval  of  25  hours. 
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plied  to  the  data.  In  retrospect,  it  is 
clear  that  the  use  of  only  one  or 
two  flying  clocks  would  have  substan¬ 
tially  decreased  the  feasibility  of  our 
experiments. 

A  hypothetical  example  may  help  to 
follow  the  logic  of  this  method.  Of  two 
clocks,  Cx  and  C>,  suppose  clock  Cx  is 
chosen  as  the  reference  clock  and  that 
at  successive  readings  of  Cx  the  time 
differences,  t.2—tx,  are  plotted  on  a 
graph.  A  change  in  slope  in  the  result¬ 
ing  time  trace  represents  a  change  in 
rate,  but  with  only  two  clocks  there 
is  no  way  of  telling  whether  C1  or  C2 
actually  changed.  But  suppose  the  time 
differences  for  a  third  clock,  are 

added  to  the  plot.  If  at  a  certain  instant 
there  is  a  correlated  change  in  slope, 
that  is,  there  occurs  in  both  traces  a 
change  in  slope  by  the  same  amount  in 
the  same  direction  at  the  same  time, 
then  by  induction  one  concludes  that 
it  is  the  clock  common  to  both  compari¬ 
sons  which  actually  changed  rate  (Ct  in 
this  hypothetical  case).  Therefore,  time 
intercomparison  data  for  an  ensemble  of 
three  or  more  clocks  permit  reconstruc¬ 
tion  of  their  rate  change  behavior  rela¬ 
tive  to  MEAN(USNO)  without  actual 
comparisons  with  MEAN(USNO),  and 
the  greater  the  number  of  clocks  in  the 
ensemble  the  more  redundant  and  self- 
consistent  the  procedure  becomes.  Rela¬ 
tivistic  effects  are  not  directly  evident  in 
the  internal  time  intercomparison  data 
because  such  changes  apply  equally  to 
all  members  of  the  flying  ensemble. 

We  recorded  the  differences  in  the 
times  indicated  by  each  member  of 
the  flying  ensemble  at  regular  intervals 
before,  during  and  after  each  trip,  that 
is,  throughout  the  entire  data  period.  An 
analysis  of  these  data  revealed  the  times 
and  magnitudes  for  correlated  rate 
changes  during  each  trip.  Thus  signif¬ 
icant  rate  changes  were  identified  and 
ascribed  to  each  clock  (9).  A  piecewise 
extrapolation  of  the  time  trace  for  each 
clock  relative  to  MEAN(USNO),  with 
proper  accounting  for  these  identified 
rate  changes,  then  produced  the  relati¬ 
vistic  time  differences  listed  in  Table  1. 

The  agreement  between  the  mean  of 
the  measured  values  and  the  predicted 
values  in  Table  1  is  very  satisfactory. 
In  addition,  the  consistency  among  the 
measured  values  is  striking.  For  the 
westward  trip,  the  standard  deviation  is 
less  than  5  percent  of  the  mean.  Ac¬ 
cording  to  the  statistical  theory  of  error, 
the  standard  deviation  is  a  valid  indica¬ 
tion  of  the  precision  of  the  measure¬ 
ment.  However,  the  number  of  mea¬ 


sured  values  is  too  small  for  a  good 
statistical  analysis. 

Another  estimate  of  the  probable 
error  in  the  final  time  differences  fol¬ 
lows  from  consideration  of  the  effects  of 
possible  errors  in  the  rates  and  rate 
changes  used  in  the  piecewise  extrap¬ 
olations.  A  preliminary  study  of  the 
effect  of  adjustments  in  the  calculated 
rates  on  the  residuals  between  the  cal¬ 
culated  and  measured  time  traces 
showed  obvious  distortions  when  de¬ 
viations  greater  than  0.4  nsec  per  hour 
were  arbitrarily  introduced.  This  value 
provides  an  estimate  of  the  maximum 
possible  error  in  the  rates  actually 
used.  The  product  of  this  rate  uncer¬ 
tainty  with  the  duration  of  the  trips 
gives  an  estimated  final  error  of  less 
than  ±30  nsec,  which  is  about  three 
times  the  standard  deviation  and  prob¬ 
ably  represents  an  overestimate. 

An  estimate  of  probable  systematic 
errors  would  require  specific  knowledge 
of  (nonrelativistic)  systematic  environ¬ 
mental  effects  on  the  clock  rates,  such 
as  a  systematic  effect  of  temperature  or 
pressure  differences.  Of  course,  any 
known  or  unknown  systematic  effect 
would  apply  equally  to  each  of  the 
flying  clocks,  for  they  all  experienced 
the  same  ambient  conditions,  so  such 
effects  would  not  be  revealed  by  the 
correlated  rate-change  analysis.  Previ¬ 
ous  studies  have  shown  that  moderate 
variations  in  temperature  and  pressure, 
such  as  those  experienced  during  the 
trips,  do  not  produce  systematic  rate 
changes.  Although  temperature  or  pres¬ 
sure  changes  sometimes  cause  individual 
random  and  unpredictable  changes  (10), 
such  random  and  uncorrelated  changes 
in  rate  do  not  contribute  to  systematic 
errors.  The  clocks  are  highly  resistant 
to  impulse  accelerations  (10).  They  are 
triply  shielded  against  the  earth’s  mag¬ 
netic  field,  and  we  found  no  effect  of 
different  orientations  in  this  field  in  the 
laboratory  (5).  In  fact  we  have  been 
unable  to  discover,  either  by  theoretical 
considerations  or  by  independent  experi¬ 
mental  tests,  any  sources  of  significant 
systematic  errors. 

In  conclusion,  we  have  shown  that 
the  effects  of  travel  on  the  time  record¬ 
ing  behavior  of  macroscopic  clocks  are 
in  reasonable  accord  with  predictions 
of  the  conventional  theory  of  relativity, 
and  that  they  can  be  observed  in  a 
straightforward  and  unambiguous  man¬ 
ner  with  relatively  inexpensive  commer¬ 
cial  jet  flights  and  commercially  avail¬ 
able  cesium  beam  clocks.  In  fact,  the 
experiments  were  so  successful  that  it 


is  not  unrealistic  to  consider  improved 
versions  designed  to  investigate  aspects 
of  the  theory  that  were  ignored  in  the 
predicted  relativistic  time  differences 
(1).  In  any  event,  there  seems  to  be 
little  basis  for  further  arguments  about 
whether  clocks  will  indicate  the  same 
time  after  a  round  trip,  for  we  find 
that  they  do  not. 
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